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We investigate experimentally an exotic state of electronic matter obtained by fine-tuning to a quantum critical 
point (QCP), realized in a spin-polarized resonant level coupled to strongly dissipative electrodes. Several 
transport scaling laws near and far from equilibrium are measured, and then accounted for theoretically. Our 
analysis reveals a splitting of the resonant level into two quasi-independent Majorana modes, one strongly 
hybridized to the leads, and the other tightly bound to the quantum dot. Residual interactions involving these 
Majorana fermions result in the observation of a striking quasi-linear non-Fermi liquid scattering rate at the 
QCP. Our devices constitute a viable alternative to topological superconductors as a platform for studying strong 
correlation effects within Majorana physics. 



Quantum phase transitions (QPT) — the singular change be- 
tween two distinct ground states driven by a non-thermal 
control parameter — are currently attracting strong interest in 
widely different fields of physics, ranging from quantum mag- 
nets and strongly correlated materials (T) to, more recently, 
cold atoms [2], nanostructures [3, 4|, and particle physics J3J. 
The foremost remarkable property of QPTs is the possibil- 
ity to create exotic quantum states of matter at the quantum 
critical point, such as deviations from the standard Fermi liq- 
uid paradigm for metals; these exotic zero temperature states 
then cause anomalous physical properties at finite tempera- 
ture JT). Another intriguing aspect of QPTs is their behav- 
ior under nonequilibrium conditions, such as either a sudden 
quench driving the system non-adiabatically through the tran- 
sition [6], or a strong perturbation provided by a large cur- 
rent density as typically realized in nanoelectronic devices 
0. Despite the ubiquity of QPTs in contemporary theoret- 
ical physics, obtaining clear experimental signatures has been 
challenging. Here, we present a thorough characterization of 
all facets of a QPT in a fully-tunable single-molecule transis- 
tor [8][9l built from a spin-polarized carbon nanotube quantum 
dot connected to strongly dissipative contacts. 

We probe the regime of resonant transport through a sin- 
gle level using both linear and non-linear conductance mea- 
surements. As a reference point, in the standard case of good 
metallic contacts, transport occurs via a resonance of finite 
width as the electronic level of the quantum dot crosses the 
Fermi energy of the leads. In contrast, the presence of dissi- 
pation drives the conductance to zero (in the limit of vanishing 
temperature), unless one tunes the center of resonant level to 
the Fermi energy and simultaneously makes the tunnel bar- 
riers between the dot and the two leads perfectly equal. In 
the latter case, the conductance saturates at the unitary limit, 
e 2 /h (e=electron charge, /i=Planck constant), as we recently 
demonstrated ||9). A quantum critical state is obtained, whose 
anomalous properties are the subject of the present study. 

By optimizing the dissipative environment, we find that the 
quantum critical properties surprisingly behave very closely 
to the result for a bound single Majorana fermion mode, caus- 




FIG. 1. Sample and schematic, a) Atomic force microscope 
(AFM) image of a sample similar to the one measured. A single car- 
bon nanotube is connected to long and resistive source (S) and drain 
(D) electrodes made of a thin Cr film. Tunnel barriers are tuned by 
two lateral side gates, b) A diagram of the quantum phase transi- 
tion in our system. The downward flow towards the quantum critical 
point (marked G = 1) is considered in Fig. [5] while the runaway 
behavior upon gate detuning is investigated in Fig. [4] 



ing in particular a quasi-linear scattering rate in temperature T 
and voltage bias V . Majorana fermions are presently the sub- 
ject of intense scrutiny in different contexts flOl . ranging from 
particle physics 1111 to topological insulators 1121 and super- 
conductors lfT3l4l5l . In our case, Majorana modes emerge due 
to two strongly interacting leads attempting to hybridize with 
the resonant level: their competition results in only partial hy- 
bridization of the level with the continuum. The quantum crit- 
ical point then corresponds to a frustrated state in which a 
non-hybridized Majorana mode is generated, as predicted in 
the related context of the two-channel Kondo model lfl6l[T7ll . 

We study the QPT by measuring the conductance along the 
flow into the QCP as well as along the runaway flow to the 
decoupled fixed point, under both equilibrium (eV <C ksT) 
and non-equilibrium (eV > fc^T) conditions. This allows us 
to present a comprehensive study of the scaling laws, extract- 
ing several critical exponents which are found to be consis- 
tently determined by a single number, the dimensionless cir- 
cuit impedance, r = Re 2 /h, reaffirming our identification of 
the quantum critical state. 

The fabrication of the carbon nanotube quantum dot was 
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FIG. 2. Differential conductance maps close to the QCP. The maps are obtained by sweeping bias V and back gate offset voltage AV^atc 
(measured relative to the center of the peak) at T = 50 mK and B = 3 T. Conductance in the symmetric case (a) is followed by maps showing 
progressively increasing barrier asymmetry in (b-d), as controlled by the side gate voltages. The scale of conductance maps in (c) and (d) 
are multiplied by x4 and x40, respectively. The prominent zero bias anomaly (ZBA) in the asymmetric cases completely disappears when 
symmetric and on resonance, indicating a QPT 



described previously J8] [9j (see Fig.[]J). The device is oper- 
ated in a magnetic field of B — 3 or 6 T and at temperatures 
down to T = 50 mK so that the quantum dot is effectively 
spin polarized. A key feature of our samples is that the cou- 
pling of the quantum dot to the left or right metallic lead can 
be varied independently using side gate voltages. We present 
data measured on two resonant levels during different mea- 
surement sessions, demonstrating compatible results. 

Differential conductance maps (G(V) = dl/dV) near one 
of these resonant peaks are shown in Fig. [2] as a function of 
bias voltage V and gate voltage offset AK at e (measured from 
the center of the peak). Throughout this paper, we express the 
differential conductance G in units of e 2 /h. Maps for differ- 
ent degrees of coupling asymmetry are shown, ranging from 
symmetric to very asymmetric coupling between the resonant 
level and the two leads. For asymmetric coupling, note the 
clear zero-bias anomaly (ZBA) — a suppressed conductance 
at zero bias for all gate values. The ZBA is eliminated only 
for symmetric tunnel couplings and only exactly on-resonance 
(Fig. [2^). This sets the location of the quantum critical point. 
These general observations [9| are summarized in Fig. [TJb) 
which sketches a quantum critical diagram for our system: at 
zero temperature, the system is in a conducting state with con- 
ductance e 2 /h at the QCP and an insulating state otherwise. A 
generic feature of QPTs JT|, which we now investigate in our 



device, is the existence of quantum critical correlations close 
to the QCP. 

We first analyze the transport data taken with symmetric 
barriers and gate voltage tuned on-resonance, which corre- 
sponds to the vertical cross section of Fig. [2^ at AV^ a te = 0. 
In this case, perfect transmission survives down to the low- 
est temperature and zero bias voltage. Fig. [3^ shows a se- 
ries of conductance vs. bias plots measured by stepping up 
the temperature. Strikingly, the conductance at the lowest 
temperatures shows a very unusual dependence on bias — a 
quasi-linear cusp. We replot 1 — G at the base temperature 
on a log-log scale in Fig. |3j). Clearly, the deviation of the 
zero-bias conductance from the unitary limit is quasi-linear: 
a fit to 1 — G oc V a yields a f=s 1.1. Essentially the same 
exponent a s» 1.2 is extracted from the T dependence of 
1 — G (Fig. BJ;). In the absence of dissipation in the leads, 
1 — G usually scales with the conventional Fermi-liquid T 2 or 
V 2 dependence. Observation of a very different, intrinsically 
non-Fermi liquid [ 16 1 exponent a underlines the dramatic role 
played by the dissipative environment. 

To account for this behavior, we draw on the analogy be- 
tween tunneling with dissipation (our situation) and tunneling 
in a Luttinger liquid - an interacting one-dimensional elec- 
tronic system ifTBI . which leads naturally to violation of the 
Fermi liquid paradigm. This analogy was formally established 
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FIG. 3. Critical flow, a) Nonlinear conductance G(V, T) at the QCP (A Vgate = and symmetric barriers) measured at increasing temperatures 
top to bottom: T = 0.055, 0.07, 0.09, 0.11, 0.15, 0.18, 0.26, 0.34, 0.41, 6.58, 0.75, 1.10, 1.36, 1.67 and 1.96 K. The quasi-linear behavior at 
low T signals the presence of a Majorana state, b) 1 — G measured vs. bias voltage V at T = 55 mK (symbols). The dashed line is a best fit to 
a power law 1 — G oc V a , excluding the range eV < ksT where the conductance saturates, c) Deviation of zero-bias conductance from the 
unitary limit 1 — G at the QCP as a function of temperature. Different symbols correspond to two different Coulomb blockade peaks measured 
on the same sample; the broken line is a best fit to Peak 2, giving an exponent a = 1.2. d) A representative set of conductance curves from 
panel (a) at four different temperatures, plotted as 1 — G(V, T) normalized by 1 — G(0, T). A line fit (see text) is superimposed on the data 
(symbols). The collapse of the data onto a single curve as a function of eV/ksT demonstrates the universality produced by proximity to the 
QCP. 



for tunneling through a single barrier in Ref. |fl8l ; we previ- 
ously extended it to the case of spinless resonant tunneling 
|9| and will further use it here (see the Supplementary Infor- 
mation). (We stress that our experiment does not probe Lut- 
tinger liquid physics in the nanotube: at this T, the length of 
an ideal, clean nanotube would have to be as large as 100 /im 
to suppress finite-size quantization.) In the bosonization lan- 
guage, commonly used to describe Luttinger liquids lfl6l . the 
approach to the QCP is controlled by the leading irrelevant op- 
erator, cos(2y / 7t6')9(? 1 19 1, which corresponds to reflection of 
electron waves from the charge accumulated at the resonant 
level. Based on the scaling dimension of this operator, and 
transcribing the Luttinger interaction parameter g to the di- 
mensionless dissipation strength r according to g = l/(r+l) 
[0Q1D, we expect the approach to the unitary limit to follow 



1 — G oc= T 2 ^ r+1 \ For this sample, the strength of the dis- 
sipation is r s» 0.75 (or g rj 0.57), as measured from the zero 
bias anomalies far off resonance (9). The predicted quantum 
critical exponent, 2/(r + 1) « 1.14, is thus very close to the 
value of a found experimentally (Fig. [3]3,c), demonstrating 
the ability of the Luttinger liquid framework to explain our 
data. 

The analogy between tunneling with dissipation and tun- 
neling in a Luttinger liquid can provide further insight into the 
microscopic nature of the quantum critical state. In the Lut- 
tinger liquid case, it is known that for g = 1/2 the resonant 
level can be mapped [20| onto the quantum critical 2-channel 
Kondo Hamiltonian |J3] |T6j |T7] . The latter may in turn be un- 
derstood in terms of a hybridization of the leads with only 
"half" of the spinless fermion residing on the resonant level, 
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leaving a fractional degeneracy associated with the remaining 
non-hybridized single Majorana mode [17|. Using this intu- 
ition, we directly map our problem (resonant tunneling with 
dissipation) for the case r = 1 and on-resonance (AV g . dte = 0) 
onto the following effective Hamiltonian: 

^Majorana = (V S - V D ) [ty(x = 0) - $f(x = 0)] Q (1) 

+ i(V s + V D )[j>){x = 0) + # f (x = 0)]b 
— 2irivFwl{x = 0)i[> c (x = 0)ab. 

Here, a and b are the two Majorana modes that describe a sin- 
gle fermion on the dot, d) = a — ib, and %/j c /f(x) are two 
extended fermionic fields that incorporate in a non-trivial way 
the electronic degrees of freedom of the leads and the electro- 
magnetic phase fluctuations in the circuit. (See the Supple- 
mentary Information.) 

One of the Majorana modes, b, is hybridized and "dis- 
solves" into the continuum. For symmetric coupling Vs = 
Yd, the first term in Eq. ([T]i is eliminated, and the other mode, 
a, is not hybridized, allowing for the existence of an inde- 
pendent Majorana mode on the dot. The density-density in- 
teraction [last term in Eq. ([T]l] does not destroy this inde- 
pendence ||2T1 l22l : however, this strong interaction (of or- 
der the Fermi energy) does govern the transport properties 
near the critical point. We show 11241 that perturbative treat- 
ment ll2TH23ll of this interaction yields an anomalous linear- 
in-temperature dependence, 1 — G oc T 1 (see Supplementary 
Information). 

This striking deviation from Fermi liquid theory follows 
from the qualitative difference in the correlations of the hy- 
bridized and independent Majorana modes: without the third 
term in Eq. ([T]), the correlations of the hybridized Majorana 
fermion decrease at long time, (b'(0)b(t)) oc l/t, while those 
of the non-hybridized mode do not decay, (cr (0)a(t)) oc 1. 
In the calculation of the conductance, the lack of decay factor 
l/t for the independent mode a reduces the exponent of the 
power-law dependence of G on frequency or temperature by 
one, reducing the Fermi liquid dependence T 2 to T 1 . We con- 
clude that the almost linear approach to the unitary conduc- 
tance observed in Fig. [3] signals the presence of a Majorana- 
like state. Our experiment thus demonstrates a route to using 
quantum critic ality to tailor exotic electronic states in nan- 
odevices. The rest of the paper aims to precisely check the 
consistency of this scenario, and to probe the non-equilibrium 
properties of this system. 

An important feature of our experiment is the possibility 
to perturb the quantum critical state by out of equilibrium ef- 
fects, namely in the presence of a finite voltage bias. We have 
already seen in Fig. [3]; that the dependences of 1 — G on T and 
V feature the same exponent a. Interestingly, the set of curves 
in Fig. [3^ associated with the low temperature, low-bias ap- 
proach to the QCP can be rescaled onto a single universal 
curve, by plotting vs. eV/ksT (see Fig.jijl). The 

solid line is a numerical derivative of the heuristic expression 



to 1 at V = 0. This scaling function describes the experimen- 
tal data quite well without any fitting parameters, as we use 
the previously determined exponent a = 2/(r + 1). This form 
of the scaling function is inspired by the universal scaling of 
g(v=o't) ' n re gi me °f vanishing conductance (away from 
resonance), which is given by the same functional expression 
with a replaced by 2r [9]. While the scaling at weak coupling 
is theoretically well established |25l . the scaling at the QCP, 
first reported here, still requires formal justification. We note 
that the similar behaviors of G for the weak coupling and 1 — G 
for strong coupling cannot be explained by duality of the re- 
flection and the transmission at the corresponding fixed points 
E6ll - indeed, the two limiting behaviors are determined by 
different operators. 

We now turn to the runaway behavior from the QCP toward 
the zero conductance state caused by a backgate offset volt- 
age AV gate , as sketched in Fig. [TJ). We focus on the shape 
of the resonant peak as a function of gate voltage and tem- 
perature or bias. In Fig. [4^, the peak shape is plotted as a 
function of the single parameter X = AV gatc /T r ^ r+ ^ 11261 . 
Clearly, the data sets measured at different temperatures col- 
lapse onto a single universal curve G(X). Importantly, this 
scaling applies not only for high values of conductance com- 
parable to e 2 /h (as captured by the resonance width at half 
height |9j), but also to relatively large values of X where 
G« 1. The scaling thereby describes the full runaway flow 
from the strongly coupled QCP to the weakly-coupled regime. 
Such single-parameter scaling has been conjectured for reso- 
nant tunneling in the Luttinger liquid in Ref. [26], where the 
double-barrier structure is treated as a single barrier with a 
gate-tunable transparency. The validity of this assumption has 
been analyzed for weak interactions [27, 28]; our results fur- 
ther corroborate the single-parameter scaling in the strongly 
interacting case (r ~ 1). 

Concentrating first on the low conductance tail obtained for 
large backgate detuning X ^> 1, we obtain G oc (AV'gatc) - ^ 
with an exponent /? ~ 3.4 (Fig. |4j;). This is very close to the 
value 2(r + 1) ~ 3.5 expected from the analogy to tunnel- 
ing in a Luttinger liquid 11261 . As in the case of the non-Fermi 
liquid exponent a found in Fig.[3]for the resonant case, this re- 
sult is markedly different from the tail of the usual Lorentzian 
lineshape, G oc (AV^ a t G )~ 2 , expected for resonant tunneling 
without dissipation. 

To further examine the assumption of single-barrier-like 
scaling, we compare the universal scaling curve obtained in 
the experiment to the recent non-perturbative renormalization 
group calculations for single-barrier tunneling in a Luttinger 
liquid |29| (Fig. [4}}). In the dissipation language, we use the 
following ansatz for the full crossover function describing the 
equilibrium conductance G as a function of both detuning and 
temperature: 



(G/G y 
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with respect to V, properly normalized with Gq = £ 2/ h and T * [AV gate ] oc ( A V g3te ) ^ . The excel- 



5 




FIG. 4. Runaway flow, a) Resonant peak conductance as a function of the one-parameter scaling variable X — AV^te/T • Several 

curves measured at eight different temperatures between 55 and 340 mK scale to the same universal curve, b) Single parameter scaling 
function extracted from panel (a) compared to the universal single-barrier scaling curve, Eq. |2](, predicted in Ref. |29|. A good fit is achieved 
by optimizing a single fitting parameter, namely a dimensionless prefactor accompanying X. (For clarity, we plot only the data measured at the 
lowest temperature.) c) Comparison of the peak shape measured at finite T under equilibrium (eV <C fesT) and non-equilibrium (eV > fc_gT) 
conditions. The horizontal axis is X = AV gal e/T r/(r+1) in the former case and X = AV ga ,e/(V/ J 4) r/(r+1) in the latter (for the value of 
constant A, see the main text). The power laws near zero conductance and near perfect transparency are the same, but a substantial difference 
between the two curves develops in the cross-over region, d) Behavior of 1 — G near the peak conductance at several temperatures. Dashed 
line corresponds to the result of Eq. which scales as ~ X 2 for small X. In each experimental curve, the change from saturation to ~ X 2 
behavior marks the range associated with the QCR Thus the quantum critical region (see Fig.[TJ is mapped out. 



lent agreement between the experimentally determined scal- 
ing function and the full runaway curve of Eq. Q (Fig. |4j)) 
vindicates the use of the single barrier expressions to describe 
the scaling of the double-barrier system away from the reso- 
nance. 

Applying a finite bias brings additional insight by allow- 
ing us to probe the relation between equilibrium and non- 
equilibrium scaling (Fig. |4j;). For the Luttinger liquid sys- 
tem, analysis near the decoupled fixed point ll26l predicts 
that the temperature scaling and bias scaling are proportional: 
b T /b v = 7r 2r [r(l + r)] 2 , where G(T = 0, V) ~ b v (eV) 2r 
at small bias and G(T, V = 0) ~ b T {k B T) 2r at low tem- 
perature. Thus, if we form a scaling curve using V in place 



of T, the two scaling curves should collapse onto each other 
in the small conductance limit if eV is divided by A = 
?r[r(l + r)] l / r k B ^ 2.8fc B for r = 0.75. The data in Fig.^ 
demonstrates that this holds for our system. However, the 
full scaling curves are clearly quite different: though they 
show, as expected [26|, the same power laws close to both the 
strong and weak coupling fixed points, the bias scaling func- 
tion makes a substantially sharper transition between the two. 
To the best of our knowledge, there is no firm theoretical pre- 
diction for the full crossover curve of non-linear conductance 
at arbitrary values of the dissipation parameter r. 

Finally, quantum critical correlations are expected [ 1 1 to ex- 
ist away from the critical point for temperatures larger than 
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a characteristic gate-dependent scale (see Fig.[T]i. Experi- 
mentally, we can define the boundary of the quantum crit- 
ical region from the behavior of 1 — G near the peak con- 
ductance. On general grounds, a dissipation-independent ex- 
ponent 1 — G oc (AKate) 2 is expected from the assump- 
tion of single-barrier scaling [26 1, resulting in a universal 
1 — G oc X 2 dependence at small X. This behavior is indeed 
observed at low enough temperatures in Fig. |4}l for X ~ 1 
meV/K 43 . However, the approach of 1 — G to is cut off 
near the QCP by the leading irrelevant operator, causing 1 — G 
to saturate at small X. (More directly, the effect of the lead- 
ing irrelevant operator is seen in the power laws of Fig. [3]) 
Thereby, the saturation delineates the quantum critical region; 
indeed, notice that the width of the saturation region gets pro- 
gressively narrower as the temperature is reduced. We may 
define the width of the quantum critical region in AV^ a t e by 
comparing the small corrections to full transparency due to 
the relevant operator oc [AVg a tc/T r / (r+1) ] 2 (Fig.H) to those 
due to the leading irrelevant operator oc T 2 ^ r+ ^ (Fig. [3}. 
The quantum critical region defined in this way has a simple 
shape: AV^ a t c oc T. We convert AV^atc to the actual energy 
position of the level Ae^ by multiplying by the experimentally 
determined "gate efficiency factor" of about 0.2. We then find 
Ae^/fesT 1 — 1, indicating that the boundary of the quantum 
critical region simply corresponds to the center of the con- 
ductance peak being shifted away from the Fermi energy by 
~ fcflT. While appearing natural, this result is not obvious a 
priori, because the width of the resonant peak is much larger 
than ksT and scales as a non-trivial power of temperature. 

In conclusion, we presented a comprehensive study of the 
quantum phase transition in a resonant level subject to finite 
dissipation, explicitly demonstrating the consistency of vari- 
ous scaling laws based on the Luttinger liquid analogy. This 
allowed us to identify the quantum critical point as closely re- 
lated to a single Majorana bound state, with a scattering rate 
quasi-linear in temperature. While our device, which requires 
fine-tuning of several gates and a carefully crafted dissipa- 
tive environment, is not likely to provide a useful platform for 
Majorana-based quantum computing OUll . it certainly allows 
exploration of fascinating Majorana physics in the presence 
of strong correlations, reaching regimes that may be harder to 
achieve in topological circuits. 
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I. EXPERIMENTAL DETAILS 

While several samples demonstrated qualitatively similar 
behavior (non-Fermi-liquid exponents close to the critical 
point), we chose to present a consistent set of data measured 
on one sample with r close enough to 1 . We show two single- 
electron peaks, labeled peaks 1 and 2 in Fig. 3c of the main 
text, measured during two cooldowns. Specifically, Fig. 2 and 
4c,d have been measured on Peak #1 at B = 3 T; Figures 3 
and 4a,b have been measured on Peak #2 at B = 6 T. The 
peaks happen to have similar parameters (as evidenced, e.g. 
by the similar width in Fig. 4a and c), so the data measured on 
the two peaks can be viewed almost interchangeably. 




v m v oiv) 



FIG. SI. a) Nonequilibrium conductance measured at AFgate = 
and at various coupling asymmetries. The curves a, b, c, and d 
correspond to on-resonance condition for the maps shown in Fig. 2 
of the main text, b) Conductance in the symmetrically coupled case, 
associated to curve a of panel a), measured at increasing gate voltage 
AV r gatc = 0, 0.7, 1.2, 1.6, 1.7, 2.0, 5.0 mV (top to bottom). 

Fig. SI represents the conductance G vs. bias voltage V 
data corresponding to the vertical cross-sections of the maps 
in Fig. 2 of the main text. Both panels include the curve mea- 
sured at the QCP (top curve) and several curves detuned from 
the QPC, either by the tunnel barrier asymmetry (panel A), 
or by shifting the level off-resonance with the gate voltage 
(panel B). Clearly, except for the special point of both sym- 
metric coupling and on-resonance condition, the conductance 
is suppressed at small bias. (Note that the zero-bias suppres- 



sion is limited at eV ~ fc^T, leading to a saturation of the 
conductance at small V .) The similar trends shown in both 
panel demonstrate the equivalent role of the gate detuning and 
tunnel barrier asymmetry in destroying the quantum critical 
state. 

In Fig. S2, we further analyze the data of Fig. 4d and present 
an empirical procedure to determine the width of the quantum 
critical region, alternative to the one discussed in the main 
text. We fit conductance G(X) with the eq. (2) of the main 
text, replacing X with an ad hoc expression (X 2 + Xq) 1 ' 2 . 
Thus introduced width, Xo, is the only fitting parameter for 
each curve. The extracted X is plotted in the inset and is 
further fit by a power law dependence on T. The resulting 
exponent of approximately 0.7 is reasonably close to the ex- 
pected 1/(1 + r). 




0.01 -| , I 
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FIG. S2. a) Select sets of data from Fig. 4d fit by the expression 
of eq. (2) of the main text, with X empirically replaced with (X 2 + 
Xq) 1 / 2 . Red dots indicate the fitted Xo for each curve. Dashed line: 
eq. (2) of the main text (without Xo). Xo is further plotted in the 
inset and fit by a power law dependence on T. 



n. DERIVATION OF THE HAMILTONIAN FOR A 
DISSIPATIVE SPINLESS QUANTUM DOT 

We provide here the microscopic theory of our experiment, 
eventually leading us, in the following sections, to the inter- 
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acting Majorana resonant level Hamiltonian given in Eq. (1) 
of the main text. 

Because of the large applied magnetic field, our experimen- 
tal setup can be first modeled by a spinless quantum dot, de- 
scribed by a fermion operator (ft, coupled to two conducting 
leads. The key feature of our device is the presence of an 
Ohmic dissipative environment, with sizeable quantum fluc- 
tuations of the voltage in the source and drain leads. The total 
Hamiltonian thus reads: 



H = H Dot + //Leads + Hj + i?Env> 



(1) 



where i/pot = e&eftd is the Hamiltonian of the dot with the 
energy level £d (tuned by the backgate voltage lg a t e ), and 
^Leads = Z)q=s,d J2k <^kc\ a c ka represents the electrons in the 
source (S) and drain (D) leads. The most important piece of 
Hamiltonian (JTJ is H T , which describes the tunneling between 
the dot and the leads with amplitudes Vs/d- 



k 



kS e 



-V D J2(cl D e lipD d+h.c), 

k 



(2) 

where the operators <ps/D induce phase fluctuations of the 
tunneling amplitude between the dot and the S/D lead. These 
phase operators are canonically conjugate to the operators 
Qs/d corresponding to charge fluctuations on the S/D junc- 
tions. We follow the standard approach to treat quantum tun- 
neling in the presence of a dissipative environment [ 1 1, which 
is valid for electrons propagating much slower than the elec- 
tromagnetic field [2J. 

It is useful to transform to variables related to the total 
charge on the dot. To that end, we introduce [ 1 1 two new 
phase operators, ip and ip, related to the phases <ps/D by 



<Ps = n s p + ip 
PD = K D lf - i> 



(3) 



where ks/d — Cs/d /(Cs + Cn) in terms of the capacitances 
of the dots to the source/drain contacts, Cqid- The phase ip 
is the variable conjugate to the fluctuations of charge on the 
dot Q c — Qs ~ Qd- Likewise, tp is the variable conjugate to 
Q = (C'sQd + C D Qs)/(C D + C s ). Assuming C s = C D , 
we have ips — (p/2 + ip and cpn = ip/2 — ip. 

The gate voltage fluctuations can be disregarded in our ex- 
periment because the capacitance of the gate is negligible, 
C g -C Cs/d (the opposite limit of a noisy gate coupled to a 
resonant level was considered in Refs. GHU). In fact, the cou- 
pling of the fluctuations of the total charge on the dot to the 
environment is negligible [1]. Thus, only the relative phase 
difference between the two leads remains HI 13, and the tun- 
neling Hamiltonian becomes 



k 



fcS"- 



; 2d + h.c. 



kD l 



'Td + h.c.1 



(4) 

The last part of Eq. ([TJ is the Hamiltonian of the environ- 
ment, i?Env (DISH- The environmental modes are repre- 
sented by harmonic oscillators described by inductances and 



capacitances such that their frequencies are given by ojk — 
l/V-^fcCfc- These oscillators are then bilinearly coupled to 
the phase operator <p through the oscillator phase: 



91 

2C 



N 

£ 

fc=i 



A. 

2Ci. 



^-{•p-'Pkf 
e J 2L k 



(5) 



III. MAPPING ONTO THE LUTTINGER LIQUID 
TUNNELING PROBLEM 



In this section, we first use bosonization [8| to map our 
model Eq. ([JJi to that of a resonant level contacted by two Lut- 
tinger liquids. In carrying out this mapping, we follow closely 
previous work on tunneling through a single barrier with an 
environment (51 [10) and the Kondo effect in the presence of 
resistive leads Q. 

The two metallic leads in our case can be reduced to two 
semi-infinite one-dimensional free fermionic baths, which are 
non-chiral [8|. By unfolding them, one can obtain two chi- 
ral fields JH, which both couple to the dot at x = 0. We 
bosonize the fermionic fields in the standard way (8] [TTI 
c S/d(x) = ^=exp[i(p s/D (x)] (neglecting Klein factors 
for simplicity), where (pgrn are bosonic fields introduced to 
describe electronic states in the leads, and a is a short time 
cutoff. Defining the flavor field <pj and charge field <f> c by 



<t>f 



">s - <Pd 



!>5 i- <PD 



we can rewrite the Hamiltonian of the leads as 



4?r 



d.v 



{d x <t>c? + l(d x <t> f ) 2 



(6) 



(7) 



with vp the Fermi velocity. The tunneling Hamiltonian then 
becomes 



V s 

H T = m eX P 

\ / 2Tra Q 



V. 



D 



\/2ira 



exp 



V2 
V2 



+ i 



d + h.c. 



d + h.c..(8) 



Note a key feature of Hj: the fields p and 0/(0) enter in the 
same way in the tunnel process. Combining these two fields 
together embodies a process which is analogous to having ef- 
fectively interacting leads as in a Luttinger liquid. 

To carry out such a recombination of the phase fields, since 
the tunneling only acts at x = 0, it is convenient to perform 
a partial trace in the partition function and integrate out fluc- 
tuations in <f) c /f(x) for all x away from x = lfT2l . For an 
Ohmic environment, one can also integrate out the harmonic 
modes at low frequency 0171 [9]. Then, the effective action 
for the leads and the environment becomes 



Leads+Env 



IT. 



\<p c (w n )\ 2 + \(f> f (u; n )f 



2r 



(9) 
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with r — Re 2 /h the dimensionless impedance of the leads 
and uj n = 2irn/(3 a Matsubara frequency. In this discrete 
representation, it is straightforward to combine the phase op- 
erator tp and the flavor field <j>f, in order to maintain canonical 
commutation relations while doing so, we use the transforma- 
tion 

1 



4>s = Vs ( <f>f + 



V2 



llr 



(10) 

(11) 



where g = 1/(1 + r) < 1. Now, the effective action for the 
leads and environment (excluding tunneling) becomes 



^Leads+Env 



K| (\Mu n )\ 2 + l^/K)| 2 + lv'K)| 2 ) 



(12) 

while the action for the tunneling part reads in the time- 
domain: 



S T = dr 



Vs 



c.c. 



D 



V / 2ttoo 



(13) 



Thus we see that the phase cp has been absorbed into the new 
flavor field <j)'^ at the expense of a modified interaction pa- 
rameter g < 1, while the new phase fluctuation cp' decouples 
from the system. It turns out that one obtains a very similar 
effective action by starting from a model of spinless resonant 
tunneling between Luttinger liquids HI 21 1131 . allowing us to 
obtain the transport properties for arbitrary r values from pre- 
vious knowledge accumulated for the Luttinger problem, as 
discussed in the main text. 



IV. QUANTUM CRITICAL POINT AS MAJORANA STATE 

In this final section, we show that the special value r = 
1, corresponding to a fine-tuned circuit impedance R = 
h/e 2 (close to our experimental value), allows a detailed 
understanding of the nature of the quantum critical point. 
Technically, a standard refermionization [8] of the tunneling 



term ( 13 i, leads to an effective interacting resonant Majorana 
level. 

The mathematical equivalence to the Majorana model starts 
by performing a unitary transformation lfl4l \l5l , U = 
exp[i(d}d — 1/2)0 C (O)/V2], to eliminate the <j) c field in the 
tunneling action, Eq. ( p"3j ), which now reads: 



St = I dr 



\J2-na^ 



d H e 

V27ra 



755*7 



+C.C. 



(14) 



This operation generates a new contact interaction between 
the dot and the phase field: 



H c = -Trv F (dU - l/2)d x <f) c {x = 0) 



(15) 



For the special value g = 1/2, we can exactly refermionize 
the problem by defining fictitious free fermion fields ip c = 
I v / 2~7rao and tpf = e^f / ^/27ra (where we neglect the 
Klein factors which do not play any role here). Electron waves 
in the leads and phase fluctuations in the circuit are thus com- 
bined into these non-interacting fermionic species. All the 
complexity of the tunneling process is now encoded in the 
following expression: 



H 



Majorana 



= Ht + H] 



Dot 



[V s 4(0)d- 



h.c. 



(16) 



[V D ip f (0)d + h. C : 



+e d dU - nvF^t(0)M°)(d^d - 1/2). 

The peculiarity of this effective Hamiltonian is the presence 
of "pairing" terms, like ipf(0)d, in contrast to the initial tun- 
neling term Eq. (|2]) where the number of the fermions is con- 
served. This structure motivates the introduction of a Majo- 
rana description of the local level, = a — ib where a and 
b are real fermions satisfying {a, b} = and a 2 = b 2 = 1/4. 
The effective Hamiltonian fTo*} is then readily expressed as 



^Majorana = (V S - V D ) ty} (0) " ^/(0)] O 

+i(V s + V D )[ip\(0)+iP f (0)]b 
+2ie d ab - 2iirv F ^t (0)^ c (0)a6; 

we have, thus, arrived at Eq. (1) of the main text. 
A special point of Hamiltonian 



(17) 



17 1 can then be identified 



when Vs = Vd (symmetric tunneling amplitudes to source 
and drain) and e d = (resonant condition), in which case 
the a Majorana does not hybridize to either the leads or the b 
Majorana level; the latter is on the other hand tunnel coupled 
to the fermion bath. In the absence of the density interac- 



tion to the field ip c (last term in Eq. ( 17 i), one ends up with 
a non-interacting Majorana resonant level model, describing 
a critical state with fractional degeneracy due to the perfect 
decoupling of the Majorana a mode (the ground state entropy 
is then S = \ log 2). Gogolin and Komnik ifTBI introduce an 
extra, fine-tuned Coulomb interaction between the quantum 
dot and the leads to get rid of the Majorana interaction term. 
However, this simplifying hypothesis is not appropriate in our 
case and leads to an incorrect description of the quantum crit- 
ical point. Indeed, for non-interacting Majoranas, the con- 
ductance at low temperature assumes a Fermi liquid behav- 
ior with quadratic in T approach to the unitary value, which 
corresponds to an exact and unfortunate cancellation of the 
leading irrelevant operator discussed in the main text. Hav- 
ing a finite Majorana interaction presents two interesting fea- 
tures: (i) the entropy of the critical state remains unchanged 
at S = s log 2 (the a mode is not hybridized with the contin- 
uum despite the density-density interaction ltT6HT8l ); (ii) an 
anomalous non-Fermi liquid scattering rate is generated, with 
linear behavior in temperature and voltage. This recovers the 
correct behavior due to the leading irrelevant operator in the 
approach to the critical point |fl9l . 



Note that the four-fermion interaction term in Eq. (17i is 
too large to be quantitatively captured by perturbation theory, 
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but rather drives the Majorana-level model into a strong cou- 
pling regime lEfl with universal scaling relations describing 
the quantum critical state in our system. Nevertheless, in- 
sight into the striking non-Fermi liquid behavior can be gained 
by calculating the conductance perturbatively in the interac- 
tion |[T9l , as mentioned in the main text. The hybridized Ma- 
jorana fermion b has a decreasing correlation function at long 
time, namely G b {t) = (b + (0)b(t)) oc l/t, identical to the 
behavior of the local (one-dimensional) bath Green's func- 
tion G^t x —Q\(t) oc l/t. However the non-hybridized mode a 
does not decay at long-times, G a (t) — (a + (0)a(t)) oc 1. The 
role of the interaction given by the last term in Eq. ( [17} is to 
generate a self-energy for the hybridized mode b; at lowest or- 
der in perturbation theory, it reads Sb(t) oc VpG a (t)G^ (t) oc 
l/t 2 . Now, the current-current correlation function of the 
original problem turns out to be proportional to the self- 
energy |19|. Its Fourier transform yields a dependence 
linear in frequency, and in this way a correction to the uni- 
tary conductance which is linear in temperature is produced, 
1 — G oc T 1 . Away from the QCP, the Majorana mode a be- 
comes hybridized, and so decays as l/t 3 , which by Fourier 
transform leads to the expected quadratic-in-energy Fermi liq- 
uid behavior. 

Finally, for r close to but not exactly 1, corresponding to 
our experimental situation (the circuit impedance cannot be 
tuned in situ but rather is fixed by the lithographic process), 



one ends up with a Majorana model similar to Eq. (17i but 



with weakly interacting Luttinger leads [ 2~TH23ll with effective 
Luttinger parameter g oc r. This residual interaction within 
the continuum leads to a slight deviation from the linear ap- 
proach to the unitary conductance, as we indeed observe in our 
data (see Fig. 2a of the main text). Although the critical state 
is not precisely a Majorana bound state in that case, the asso- 
ciated ground state still possesses entropy S — \ log(l + r) 
associated with a non-trivial fractional degeneracy lfT71l23l . 
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